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We have recently developed a canonical transformation of the Hubbard 
and related models, valid for systems of arbitrary size and for the full 
plane; this is particularly suited to study hole pairing. In this work we 
show that exact diagonalization results of the one band Hubbard model 
for small clusters with periodic boundary conditions agree well with the 
analytical ones obtained by means of our canonical transformation. In 
the presence of a pairing instability, the analytic approach allows us to 
identify the Cooper pairs. They are W=0 pairs, that is, singlet two-hole 
eigenstates of the Hubbard Hamiltonian with vanishing on-site repulsion. 
Indeed, we find that the Coulomb interaction effects on W=0 pairs are 
dynamically small, and repulsive or attractive, depending on the filling. 



1 Introduction: the Hubbard model 

Let us consider the Hubbard model on a square lattice of N x N sites, with 
hamiltonian 
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H = H + W, 



(1) 



where 




(2) 



<T (r,r') 



with r, r' nearest neighbours and 



W = U^n r] n rl . (3) 

r 

We use periodic boundary conditions and the one-body wave vectors are k = 
(k x , k y ) — tJKp, q) with p and q integers. 

2 W=0 pairs 

The strong on-site repulsion between two opposite spin fermions normally pre- 
vents the formation of singlet bound states. However, the planar square sym- 
metry of the Hubbard and related models allows singlet two body eigenstates of 
H belonging to the zero eigenvalue of W. We will refer to these states as W=0 
pairs. This means that the two fermions of a W=0 pair do not interact directly, 
but only by means of virtual electron-hole excitations. The background parti- 
cles play a crucial role in determining the structure of the effective interaction 
which can be, in priciple, attractive or repulsive. 

In previous works 1 ' 2 we have shown how to build W=0 pairs with zero total 
momentum. There the idea was to project the determinantal state 

\d(k)} = c^d kl \vac) (4) 

on the Irreducible Representations (irreps) of the point symmetry group C± v . 
Remarkably, by projecting on the irreps A 2l B\ and B 2 one obtains exclusively 
W=0 pairs. Here we want to point out a new and more general criterion to get 
W=0 pairs. Let Q a symmetry group of the non interacting Hubbard Hamil- 
tonian Ho, big enough to justify the degeneracy of the single particle energy 
levels. Let us consider a two body state of opposite spins trasforming as the 
i-th component of the irrep L of Q: 

I^HE^i'^c^Jmc) (5) 

Then we have 

n r ^n rl \^P) = *P(r,r)4 T 4j«ac) = *P(r,r)\r |,r |). (6) 

Let be the projection operator on the i-th component of the irrep T. Since 
^(r) £ ^(r) (r , r) |r T ) r |} = ^ ( r , r) | r T , r |) ( 7 ) 

r r 

if 

ffVT,r|} = Vr (8) 

then 

*f } (r,r)=0 Vr. (9) 



Clearly eq.(||) is true if and only if 

PP\ra)=0 Vr 

where \ra) = c\ a \vac). 

It is always possible to write \ra) as 

where £ is the set of the irreps of the one-body spectrum of Hq and \ip\ g .) the 
corresponding eigenstate with spin a. From ([h]) it follows directly that if V 
does not belong to £ 

P ir,) \ra) = 

and so 

P {F ">\rtr I) =0. 

We have proven the following 

THEOREM: Let |W) be a two-body eigenstate of the kinetic energy 
Hq with spin S z — 0. Projecting on an irrep not contained in £ , 
we get either zero or an eigenstate of Hq with no double occupancy. 

The singlet component of this state is a W=0 pair, and its momentum does 
not generally vanish. 

3 Binding energy: analytic approach 

From now on we call hole the fermion created by and electron its antifermion. 

The Schrodinger equation for the ground state of our N x N square lattice 
with iih holes is 

H\9 ) =E h (n h )\^ )- (11) 

If the non interacting ground state with — 2 holes can be written in terms of 
a single determinantal state, the exact |^o) can always be expanded in terms 
of excitations over it: 

l*o) = E a ™i m ) +E & «' a ) + ^2 C P\0) + — ( l2 ) 

m a {3 

here m runs over pair states, a over 4-body states (2 holes and 1 electron-hole 
pair), (3 over 6-body ones (2 holes and 2 electron- hole pairs) and so on. Eq.([l2"|) 
is an expansion in the number of virtual excitations. 

In the following we set up a procedure which 1) carries out the "Configura- 
tion Interaction" calculation in a compact and efficient way 2) separates neatly 
the effective interaction from the self-energy contributions to the ground state 
energy. To understand how this mechanism works, for the moment we truncate 
the expansion to the (3 states. Then equation ([ll]) yields 

(E m - E h (n h )) a m + ^a m >W m . m > + }b a W miCt + )cpW m ,p = (13) 



(E a - E h (n h )) b a + Y^ a m' W a , m > + Yy°<' W »^ + ^cpW a ,p = (14) 

m' a* {3 

(Ep - E h (n h )) c + ^a m /W^, ro / + ^bcfWp,*' + ^c^W^ = 0. (15) 

m' a.' j3' 

Choosing the /3 states in such a way that 

(Ho + W) M , = E'p8 m 
we exactly decouple the 6-body states getting 

(E m - E h (n h )) a m + ^a m W m m , + ^6 Q W^ Q = 

m' a. 

(E a - E h (n h )) b a + 5^0m'W£ )m , + Yy^' W '^ = 
m' a' 

where W"s are the renormalized interaction coefficients. It is clear that if wc 
had truncated the expansion to an arbitrary number n of electron-hole virtual 
exitations, we should have obtained the same results but with further renor- 
malizations. This is a recursion method to perform the full canonical transfor- 
mation; it applies to all the higher order interactions, and we can recast our 
problem as if only 2 and 4-body states existed. Now we choose the a states in 
such a way that 

(H + W') a ,a, = E'J a , a , 
getting the following eigenvalue problem 

(E h (n h )-E m )a m =J2 a m' (m\F[E h (n h )] + W eff [E h (n h )]\m') , (16) 

m 1 

where 

(m\F[E h (n h )]+W eff [E h (n h )]\m') = W' m , m , + £ H^HH^V') _ 

Equation ( [l6| ) determines the amplitudes a m of the m states in the nh hole 
ground state and the corresponding eigenvalue Eh(rih) relative to the hole vac- 
uum. Here F is the forward scattering operator and W e ff the effective inter- 
action. After the shift 

Ho^H -E h (n h -2), (18) 

it is perfectly consistent to interpret a m as the wave function of the dressed 
pair. It is worth to underline that this canonical transformation enables us to 
identify the effective interaction between two holes ina%-2 holes background 
self-consistently. In principle we can work it out analytically and the expansion 
is neither in U nor in t but in the number of vitual exitations. It permits us 
to identify the binding energy |A(n/j)| too. After the shift ( |l8| ) we write the 
ground state energy Eh(nh) of the Hubbard hamiltonian as 



E h (n h ) = 2E F + A(n h ) 



(19) 



where Ef is the renormalized Fermi energy and A(n^) the two holes energy 
gain of the interacting system with respect to the non interacting one. If the 
effective interaction W e ff is attractive and A(n^) < we can speak of hole 
pairing. 

4 Pairing in the 4x4 lattice Hubbard model 

In this section we want to compare the results for the binding energy |A(n^)| 
obtained 1) from eq.([l9|) by truncating the canonical transformation to the a 
states and 2) from the definition 

A{n h ) = E h {n h ) + E h (n h - 2) - 2E h {n h - 1) (20) 

by computing exactly the three ground state energies involved in eq.(p0|). The 
equivalence of the two definitions of the binding energy was shown in Refs [2] , [6] . 

We consider a 4x4 lattice Hubbard model with t = — 1 and periodic bound- 
ary conditions; we investigate the possibility of pairing when two holes are added 
to a two hole background. The one particle energy spectrum of Hq has 5 equally 
spaced levels having degeneracy 1,4,6,4,1 respectively. 

Our first task is to determine of the symmetry group Q. The Space Group 
containing the translations and the 8 C± v operations is not enough to explain 
the degeneracy 6. Indeed the largest dimension of the irreps of the Space Group 
for this 4x4 lattice is 4. The Space Group has 128 elements and 20 classes. 
As observed by previous authors 3 ' 4 there must be additional space symme- 
tries. We have found how to obtain them. Label the sites from 1 to 16 in 
the natural way; then rotate the plaquettes 1,2,5,6 and 11,12,15,16 clockwise 
and the other two counterclockwise by 90 degrees. This is one of the missing 
space symmetries, as one can see that each site preserves its first neighbours. 
It is not an isometry. The other missing symmetries are generated by adding 
this one to the Space Group. In this way we obtain the Symmetry Group Q 
of 384 elements, which is large enough to justify the degeneracy of each one 
particle energy level. It still has 20 classes; the dimensions of the irreps are 
1,1,1,1,2,2,3,3,3,3,4,4,4,4,6,6,6,6,8,8. 

The exact interacting ground state of the system with 4 holes is threefold 
degenerate and belongs to an irrep of Q that contains the irreps Ai and Bi of 
Ci v once and two times respectively. This irrep is not contained in the one 
particle spectrum of Hq, which does not admit degeneracy 3. By the above 
theorem, pairs belonging to it must be W=0 pairs. These W=0 pairs arise from 
the single-particle states of Hq with eigenvalue -2, and their irrep is contained 
in the square of the irrep of the single-particle level. 

The m states of the expansion (|lj) are all the W=0 pairs belonging to the 
ground state irrep. We have computed the binding energy |A(4)| analitically 
from eq. ( ft9| ) by truncating the expansion (|l2| ) to the a states and also by means 
of exact diagonalization from eq. (^0|) . The results are listed below for different 
values of the on-site interaction U (energies are in eV): 





U =0.1 


U =0.5 


U =0.7 


U =1 


^(4) exact 


-0.053 


-0.32 


-0.67 


-1.18 


^(4) analytic 


-0.078 


-1.95 


-3.83 


-7.81 



As expected, with increasing U the difference |A(4) ea;act — ana i yt ic\ increases 
because the renormalization induced by virtual electron-hole exitations becomes 
important and it is no longer a good approximation to consider the a states only. 
Nevertheless the canonical transformation still predicts the right sign of A. 

The above canonical transformation applies when two holes are added to a 
determinantal vacuum. To study the system at and close to half filling, we have 
extended the above canonical transformation to the case when the vacuum is de- 
generate. We are able to demonstrate analytically that holes pair and to explain 
the ground state symmetries that were found numerically in Rcf.[3],[5]. Thus, 
W=0 pairs are responsible for pairing and for the symmetry of the interacting 
ground state. This is also confirmed by the superconducting quantization of 
magnetic flux; in previous works 2 ' 6 we show that the C± v symmetry is restored 
exactly at half fluxon 0o/2 = hc/2e allowing the existence of W=0 pairs. The 
symmetry of the interacting ground state is still the same of the possible W=0 
pairs and the corresponding energy versus flux has there a second minimum. 
A full account of the theory will be submitted elsewhere 7 . It is clear by now, 
however, that the existence of bound pairs of nonvanihing momentum opens up 
the possibility of a Jahn- Teller distortion and the present approach is likely to 
predict charge inomogeneities. 
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